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ABSTRACT 

Chromospheric wave activity around flares and filaments has been a research 
focus for years, and could provide indirect measurements of local conditions that 
are not otherwise accessible. One interesting observed phenomenon is oscillations 
in filaments, activated by distant flares and the large-scale waves they produce. 
Characteristics of these oscillations, such as periods, amplitudes, and lifetimes, 
can provide unique information about the filament. 

We measure oscillation properties in flares and filaments from Ha chromo- 
spheric data using a new method that provides important spatial and frequency 
content of the dynamics. We apply the method to two flare events where filaments 
are observed to oscillate and determine their properties. 

We find strong oscillatory signal in flaring active regions in the chromosphere 
over a range of frequencies. Two filaments are found to oscillate without any 
detectable chromospheric wave acting as an activation mechanism. We find that 
filaments oscillate with periods of tens of minutes, but variations are significant 
at small spatial scales along the filamentary region. 

The results suggest there is a frequency dependence of the oscillation ampli- 
tude, as well as a spatial dependence along single filaments that is more difficult 
to quantify. It also appears that the strength of the oscillations do not necessarily 
depend on the strength of the trigger, although there are other possible effects 
that make this conclusion preliminary. Applications of this technique to other 
events and different data sets will provide important new insights into the lo- 
cal energy densities and magnetic fields associated with dynamic chromospheric 
structures. 



Subject headings: Sun: atmosphere - Sun: oscillations - Sun: filaments, promi- 
nences - Sun: flares - techniques: image processing 
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Introduction 



It has been known for a long time (IRamsey fc Smith! Il966l ) that oscillations in solar 
prominences and filaments can be excited by distant flares. Such events have generally been 
characterized by so-called large- and small-amplitude oscillations according to the amplitude 
of the their observed velocities (compared to the chromospheric sound and Alfven velocities). 
Large- amplitude filament oscillations display speeds upwards of 20 km s -1 and periods of tens 
of minutes to hours, while small-amplitude oscillations typ ically are a few kms " 1 , with a less 
well-constrained range of periods, from minutes to hours (lArregui et al.ll2012l ). Motions are 
found in both transverse and longitudinal directions with respect to the filament axis. 

The chromosphere displays a wealth of wav e phenomena. Earl y detections o f oscil- 
lations in Ha measurements were carried out by lElliottl (119691 ). and lHarvey et al.l (119931 ) 
used several chromospheric lines to observe the high-frequency (5mHz) p modes, as well 
as 5-min. acoustic modes, both escaping from the photosphere. However, a flare (or some 
other source) activating a filament oscillation is a rather rare phenomenon, quite different 
than the oscillations that are always present in the chromosphere. Triggers of the filament 
activation are thought to be Moreton waves or EUV waves produced from remote flares 



(EtoetaL 



2002 



Okamoto et al. 2004; Balasubramaniam et a 



Li &: Zhand I2012T). or iets and su bflares (jVrsnak et al.l 120071 ). or even erupting filaments 



2007. 2010; Asai et al. 2012 



themselves (llsobe fc Tripathill2006l ). It is likely that intrinsic properties of the chromosphere 
can be deduced from studying these oscillations, as filaments seem to have associated eigen- 
frequencies when excited. This could be used to probe their magnetic structure and topology 
more precisely than traditional observations. 



Recent reviews (ITripathi et al.ll2009l ; lArregui et al.ll2012l ) point out that much of the de- 
tailed dynamics between the interaction of the Moreton or EUV waves with the filament, the 
oscillatory structure (vertical/horizontal and longitudinal/transverse), the damping mecha- 
nism, and possible causes of filament eruptions are still very open issues, as these filament 
oscillations have only been observed in about a dozen cases. 

The scope of this work is not to directly address all the issues above, but to provide a new 
tool for analyzing these and similar oscillatory phenomena. We have ob served two osc i llatin g 
filament events, one never detected and the other previously studied by lGilbert et al.l ( 120081 ). 
Standard analysis in previous studies typically involves reducing a three-dimensional (3D) 
data cube to a one-dimensional (ID) time series by spatially averaging over pixels in the 
filament that appear to be oscillating, and then performing a wavelet transform in time- 
frequency space to ascertain the variable frequency content. The reduction in dimensions 
has the effect of smearing possibly important spatial information. At the same time, the 
result may be strongly dependent on how the averaging is carried out. 
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Here we describe a new technique for analyzing phenomena that have dynamic spatial 
and temporal scales. It is an extension of a wavelet-type technique, preserving the 3D nature 
of the input data. Upon application, we observe properties of chromospheric oscillations that 
do not easily fit into the two categories described above. We observe important differences in 
flaring regions and regions of oscillating filaments, indicating different responses to the trigger 
mechanism. We also find anomalies in oscillating filaments that have and have not interacted 
with a strong Moreton wave. A standard analysis is also carried out for comparison. 

To introduce this technique, we first describe the example data sets in § [21 We provide 
an overview of the algorithm in § O which is then applied to the data and results are shown 
in § HI After comparisons of the results to standard methods and a summary are provided, 
discussion is given in § [5J 



Data 



We study two flare events observed with the United States Air Force/National So- 
lar Observatory Improv ed Solar Observing Optical Network (ISOON) telescope located at 
Sacramento Peak, NM (INeidig et al.lll998l ). Full-disk solar images (2048 x 2048 pixels at 
1.1" sampling) at 1-minute cadence are observed in the 6563 A Ha line center with a 0.1 A 
bandpass filter. Corresponding Doppler maps are also utilized, and have been constructed 
by the following procedure: 14 images across 14 wavelengths centered around the spectral 
line were gathered on a day in 2004 when there was little activity on the Sun. After cali- 
brating the images for intensity corrections and flat-fields, spectral line profiles for each 1.1 
arcsecond pixel were obtained, separately. From the center-of-gravity (spectral line center) 
for each spectral line, the intensity-difference at ±0.4A, and a Doppler shift (line-center for 
each pixel, mean line-center for entire image) was measured. By plotting all the points of 
Doppler shift versus the intensity difference, a characteristic relationship is found. For any 
day of observing the intensity differences, a conversion of that measure to Doppler shift is 
computed (assuming that the calibration hol ds true). In the flares and filam ents considered 
here the Ha line wings remain in absorption (jBalasubramaniam et al.l 120041 ) . 



A solar flare of x-ray strength M8.0 occurred on May 13, 2005 around 16:13UT in active 



region NOAA 10759. It was a classic two-ribbon flare that also cor respondec 



CME that delivered the most intense geomagnetic activity of 2005 (IBisi et al. 



characteristics of this well-kno wn event have been studied (jYurchyshyn et al 



to a very large 



20JL0I )- Various 



20061: Liu et al. 



20071 : iKazachenko et al.ll2009l ). although this is the first time oscillations in the active region 
and the nearby filaments have been detected. The second event occured on December 6, 
2006, where an X-class flare erupted at about 18:40UT in NOAA 10930. For both cases, 
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ISOON provided high-quality intensity and Doppler data from about 4 hours before until 
4 hours after flare initiation. Pre-flare snapshots of full-disk Ha intensity images for each 
data set are shown in Fig. [TJ The flaring active regions (bright in Ha), as well as several 
quiescent filaments (dark in Ha) that this investigation targets are outlined in boxes. 

One important differ ence between these two flares is that the 2006 ey ent produced a 
prominent Moreton wave ([Gilbert et al.l 120081 ; iBalasubramaniam et al.l 120101 ) , a fast-moving 
chromospheric disturbance. This particular wave propagated from the active region and 
interacted with several filaments on the solar disk, some quite far fro m the excitation loca tion. 
Filament- wave interactions in this flare were studied in detail by iGilbert et al.l ( 120081 ). A 
second important difference is the active region responsible for the 2006 event is located close 
to the limb, where strong (horizontal) velocity signals should be apparent in the line-of-sight 
Doppler observations. 



3. Data Analysis 

Studies of oscillations in filaments in the past have used velocity time series maps from 
both wings of a spectral line and determine when the filament transitions from red- to 
blueshifted, giving an estimation of the period (for example, Okamoto et al.|[2004 ; Gilbert et al 
20081 ) . Frequency content has also been derived by spatially averaging a set of pixels in the 
region of interest to produce a proxy ID "light curve," and then usi ng wavelet techniques 
J Jain k Tripathv!ll998l : iMaurva k Ambasthall2008l ; IPinter et al1l2008h . 



We demonstrate a simple and powerful way of visualizing these types of data, that pre- 
serves both frequency and spatio-temporal information. We generate what we call " frequency- 
filtered amplitude movies" (FFAMs), with a method similar to what would be described as 
3D wavelet analysis. 

To produce the FFAMs, consider a general data cube (we will be working with intensity 
and Doppler velocity) of 2 spatial and 1 temporal dimensions, f(x,y,t). Assume there are 
N time steps (images) available. For each time step ti, the following steps are carried out: 



1. Filter the data in the time domain over a segment of length T, which in this study 
is always taken to be 60 minutes, thus comprising Nt = 60 images. The 60-min 
filter tapers off quickly, but smoothly to zero at each "edge." In other words, only 
the subsequent Nt images after the current time step are considered. The filtering is 
represented as 

fifay^T) = F^t; [U,U + T}) ■ f(x,y,t), (1) 
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where the brackets in the second argument of the filter F\ denote the range where it is 
nonzero and equal to unity. The filtered data f i are labeled by the particular starting 
time step i and the segment length T. 

2. Compute the temporal power spectrum of that segment 



where v is the cyclic frequency. Only the data over the segment T will contribute to 
the power. We subsequently convert this power spectrum to an amplitude spectrum 
Ai(x, y, v\ T) to get units of velocity when using Doppler velocity data. 

3. Apply a second filter in frequency to the amplitude spectrum, which has a central 
frequency u and a bandwidth of Au, and compute the average over the frequencies in 
that particular band: 



where the (. . .)„ operator indicates an average over the values of frequency in the range 
u ± Av/2. The amplitude quantity A, is a two-dimensional image labeled by central 
frequency, bandwidth, and segment length. 

4. Continue for all % < N — Nt (the final T minutes are not computed to avoid the 
"edge"). The final data product, FFAM = ^Aj, is a three-dimensional quantity in 
x,y,t. It can be visualized as a movie (see online examples). 

We consider frequencies in 1 mHz-wide bands, centered on the values u = 1, 3, 5, and 
7 niHz. The temporal Nyquist frequency for the ISOON cadence is 8.3 mHz. The data 
sets in this study are approximately 420 minutes in length. Since segments from which 
we compute the amplitude spectrum are only 60 min in length, the frequency sampling is 
du ~ 0.28 mHz. Therefore, the average over each 1 mHz band uses about 3-4 images. We 
have done experiments with shorter and longer T, as well as with different bandwidths and 
central frequencies, and the results are qualitatively similar to those shown in § HJ 




2 




Ai(x, y\ i/ , Au, T) = ( F 2 (u) u ± Ai//2) A^x, y, v\ T) ) 



(3) 



4. 



Results 



We compute and study FFAMs of the oscillatory dynamics associated with the ISOON 
flaring events at different frequencies, and discuss the active region oscillations and the 
filament oscillations separately. 



-6 - 



4.1. Oscillations in flare regions 

The active regions in which the flares occur in the 2005 and 2006 events show oscillatory 
features over the entire frequency range studied here. The likelihood of a strong contribution 
from the underlying photospheric acoustic p modes can be ruled out for several reasons. The 
dramatic increase in amplitude occurs only once the flare erupts. If the p modes were 
responsible, their effects would be present over the entire data sequence. Secondly, the 
frequency range considered here is mostly below the photopsheric acoustic-cutoff frequency 
(~ 5.3 niHz). The fil ter bandpass o f ISO ON is also quite narrow, much more so than 



the instrument used by lHarvey et al.l ( 119931 ). where both 5- and 3-minutes oscillations were 
seen in Ha. We are confident that the effects studied here are confined mostly to the 
chromosphere. 

Figure [2] shows oscillations as a function of frequency in one snapshot of the FFAM from 
velocity data for flare regions FR4 and FR 2 . The central time of the 60-min segment used in 
the computation is provided. The figures represents the slice at the maximal amplitude of the 
time series. For FRi, the UT 16:22 moment is about 9 minutes after the flare erupted. Thus, 
there is both pre- and post-flare signal in the derived amplitude. For FR2, the 18:17 slice is 
about 23 minutes before the flare erupted, thus there is mostly pre-flare signal involved in 
the computation. 

It is evident from the figure that the flares "excite" all frequencies from about 0.5 to 
7.5 mHz to detectable amplitudes in velocity (note the figure is on a logarithmic scale). We 
are unaware of any previous studies that demonstrate this. The velocities in the 1 mHz band 
reach about 3kms _1 and 4 km s" 1 at peak amplitude for 2005 and 2006 data, respectively, 
and are usually reduced at higher frequencies. The stronger amplitudes in FR 2 could be due 
to the fact that it is near the limb, thus implying stronger horizontal motions. Interestingly, 
the 2006 flare's maximal amplitude is in the 3 mHz band, at about 5kms _1 . In intensity 
the higher frequencies near each flare contribute much less oscillatory power, as is discussed 
further in Sec. 14.31 

Once the flare erupts, the excess oscillatory amplitude lasts for about 100 minutes after 
the peak in the FFAM, although this strongly depends on the location near the ribbons. 
Some regions are damped much faster or more strongly than others. Due to the 60 min 
time-averaging involved in creating the FFAMs, it is not convenient to determine the precise 
duration of the increase in oscillation amplitude, as the time resolution is insufficient to give 
certainty to this number. More spatial structure is seen in the 1 mHz band, especially in the 
plage area of the active regions, and the two flare ribbons in FRi are somewhat visible in the 
oscillation amplitudes, as is the sunspot on the northern edge. FFAMs were computed using 
the ISOON Ha line core (intensity) data as well, and the same general spatial characteristics 
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are seen, except for weaker power at high frequencies. 



4.2. Low-frequency oscillations in filaments 

Throughout these two data sets the solar disk was littered with quiescent filaments. We 
identify and analyze two of them for each event, denoted by FLi_4. The FFAMs show a very 
dramatic "lighting up" (a rapid increase in oscillation amplitude) of the filaments well after 
the distant flares erupt, although each event has somewhat different characteristics. 

Figure [3] shows snapshots from the FFAM of filaments FLi and FL 2 in the 2005 time 
series. The images are those where the filament oscillation amplitude peaks in the FFAM, 
corresponding to the time segment centered at 55 and 62 minutes after the flare erupts for 
FLi and FL2, respectively. One explanation of the different timings is that the filaments are 
at different distances from the flare center (see Fig. [I]), and whatever activation mechanism 
is at work (such as a chromospheric wave) needed to travel different distances. 

The noticeable difference compared to the oscillations in the flare region, is that the 
filament only shows strong amplitude excess in the lowest frequency bandpass of 1±0.5 mHz, 
even on the logarithmic scale. This corresponds to periods in the range from about 11 to 33 
minutes. The amplitude is significantly above the background level in this passband, lasting 
for about 100 minutes. The velocities of the oscillations are on the order of 0.5 — 1.0 km s -1 . 
This is the first detection of oscillatory behavior in these two particular filaments. 



As in previous studies (ITripathi et al.ll2009l ). one might conclude that anywhere from 



about 3 to 10 periods of the filament oscillation took place after it was excited. With our 
method, it may be too simplistic to make such a conclusion, since there are small-scale spatial 
structures that are operating on different time scales, as is clearly evident in the FFAMs. 
Some regions show low-frequency power for the whole duration (like the central part of the 
filament), while others do not (the "ends" of the filament). 

The 2006 quiescent filaments are shown in Fig. HI We note again that a strong More- 
ton wave was associated with this flare that propagated ou twardly from the active region 



at speeds up to 1600 km s 1 (IBalasubramaniam et al.ll2010l ). The wave passed near both 



filaments and presumably activated their oscillations. The maximum values of the ampli- 
tudes in this case are ~ 0.2 km s" 1 , less than half the strength than found in the FLx and 
FL2 filaments. Filament FL3 shows no oscillations in its center, only on its ends, which is 
a curious feature. Additionally, even though these two filaments are not equidistant to the 
flare region, the maximal oscillation amplitude occurred at approximately the same time. 
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That the velocity oscillation signal in the filaments is weaker than in the 2005 event is 
surprising, given that there is a clear and strong Moreton wave observed that interacts with 
the filaments. Possible explanations for this discrepancy might include viewing geometry, 
intrinsic difference in filament properties, lower quality of data, or the overall distance of the 
filaments from the flare, which is larger than in the 2005 case. 



4.3. Comparison of method to other analysis techniques 

Analyzing one- dimensional time series of various regions by wavelet an alysis is one 



techn ique that has become standard for the study of filament oscillations (e.g., iPinter et al. 



20081 ). Computing "stac k plots" of the sign al along a certain slice of the data as a function 



of time is another (e.g., iLi fc Zhangl 120121 ) . We carry out such analyses to compare and 
contrast the results with those of the FFAMs for the flare FRi and filament FLi regions. 
The other regions show similar behavior. 

Figure [5] shows an analysis of FRi. A set of pixels near the flare ribbons was averaged 
and time series of intensity and velocity were obtained. The intensity light curve shows 
the characteristic flare impulse and exponential decay, accompanied by low-level oscilla- 
tions, while the velocity shows a fluctuating signal on top of a sudden blueshift (towards 
observer) as the flare erupts. Other pixels in the flare-ribbon region show redshifts at flare 
onset (for a detailed discus sion of velocity profiles in broadband Ha measurements, see 



Ichimoto fc Kurokawal Il984l ). In both cases the frequency content of the oscillations are 
difficult to quantify in the time series alone. Wavelet analysis of this set of pixels confirms 
that all frequencies in velocity are activated, but that the low frequencies are strongest, as 
the FFAMs indicated. The strongest amplitudes last for about 80 minutes. The spectrum 
in intensity does not show much high-frequency power. These conclusions change somewhat 
as different regions are analyzed. 

Perpendicular slices through the flare ribbons passing through the center of the same 
set of pixels are used to produce stack plots for intensity and velocity and are given in Fig. El 
From the horizontal slice one sees that the east side of the flare (closest to point C) is the 
first to begin oscillating, as is evident also in the FFAM. The strongest signal comes from 
the outer edge of the ribbons. The velocity stack plots indicate that the east-west direction 
(C-D) has a higher signal-to-noise, even though the other direction has the larger overall 
amplitude. 

Figure [7| shows a similar analysis for the filament FLi. A set of pixels near the center 
of the filament that showed strong oscillations in velocity were chosen and averaged to give 
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the time series. There are 6-7 long-period oscillations (~ 20 min) at amplitudes of about 
1 kms -1 in this region, and these show up clearly in the wavelet spectra. In intensity, the 
variations in the time series are less well defined, but the wavelet spectra do reveal the low- 
frequency nature of the oscillations. Contrasting with the chosen flaring region, the strong 
amplitudes near the filament persist for well over 100 minutes. The stack plots in Fig. [8] 
show slices through the filament along perpendicular directions. Again, the global view of 
the intensity slice is not particularly useful for studying any variations, while the velocity 
plots clearly show the filament being triggered by the flare and its subsequent oscillations. 
The motion appears to be primarily in the east-west longitudinal direction. 

The use of these standard techniques are indeed very useful for studying these phenom- 
ena, especially when coupled with an analysis using FFAMs. However, the choice of pixels 
and slice coordinates to produce Figs. (MS requires a lot of trial and error, and the inferences 
vary strongly from region to region. This could be expedited by first using the FFAMs to 
choose suitable regions that show interesting behavior, and then looking into more detail 
with the other methods. It is interesting that the stack plots in intensity for the filament do 
not strongly show the flare eruption, at least along these particular slices, while the intensity 
FFAMs and wavelets clearly do. 



4.4. Summary of results 



To summarize the analyses, we observe the fascinating occurence of rapid communication 
at large distances across the solar disk. In the 2005 event, strong amplitude increases on the 
far western solar limb are observed, very far from the flare area in the eastern hemisphere. 
Indeed, Fig. M shows full-disk snapshots of velocity FFAMs in the two lowest bandpasses near 
the peak flaring times. The filamentary amplitudes are evident everywhere at low frequency, 
and rather absent at 3 mHz and higher. For the 2006 flare, the 1 mHz band shows weak 
excess acoustic power just outside the flare region, and one barely discerns the brightenings 
in filaments FL 3 and FL 4 in this global picture. 



A com parison with o t her la rge-amplitude filament oscillations as reviewed by lTripathi et al 



(120091 ) and lArregui et al.l (120121 ) is useful. Other p ublish ed studies find oscillations with min 



imum amplitudes near 20 kms 1 . iGilbert et al.l (120081 ) find in line-of-sight velocity images 



from the Polarimeter for Inner Coronal Studies (PICS) Ha instrument that the FL 3 fila- 
ment oscillates at amplitudes up to 41 kms" 1 . We find amplitudes of a few kms -1 in the 
flare regions themselves, but less than lkms -1 in the filaments at low frequency. In the 
raw Doppler signal from ISOON we find maximum velocities of the FL3 filament of several 
kms -1 . 
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While we do observe features that approach several tens of kms -1 in the line-of-sight 
velocity data, particularly in the flare (see Fig. [61 for example), the coherent variations over 
an extended time period do not possess such amplitudes. We attribute the discrepancy with 
other values in the literature to the fact that the FFAMs select out the persistent variations 
over longer time scales (in the computations here, at least 60 min), and thus it is a matter 
of effective time resolution. Differences in instrumentation likely have some effects too. 

For filament osci llation periods, studies have found a range from about 6 to 150 minutes 



(ITripathi et all 120091 ). The periods found here are on the lower end of those values, but 
still consistent with them. One might argue that small-amplitude oscillations are a suitable 
classification for these events, which tend to have a few kms -1 amplitudes, but can also have 
very short periods on the order of seconds to a few minutes, which we are unable to observe. 



Gilbert et al.l (120081 ) have already studied both FL 3 and FL 4 using different data and 
standard methods. For FL3, they find 3 hr of oscillatory motions, somewhat longer than 
found here by about 30 minutes. They also find periods from 18-39 min. For FL4 they find 
the duration of oscillations to be only about 17 minutes, whereas we observe increased power 
(at 1 mHz) for more than 1 hr. The different diagnostics as a function of frequency is an 
advantage of the method presented here. 



5. Discussion 

This work is important for two reasons: (1) we have presented a simple algorithm 
for both visualizing and studying dynamic phenomena in 2+1 dimensions with differing 
frequency content; (2) we have used this technique to analyze two different events that 
exhibit low-frequency filament oscillations excited by a distant flare. We observe robust 
oscillations in the active regions where the flares erupt over a broad range of frequencies. 
Photospheric modes can be ruled out as the source since the oscillation amplitude only 
increases after the flare initiates. 

Strangely, the weaker filament oscillations are observed in the system that interacted 
with a strong Moreton wave, which can trigger of such oscillations, but it is possible that 
line-of-sight effects or wave propagation direction can be involved in these differences. The 
velocities found in the FFAM analysis do not fit into the large- and small-amplitude oscil- 
lation classifications. It might be useful to have sub-categories of such events that show 
oscillating strengths as a function of frequency. We have been able to confirm the work of 



Gilbert et al.l ( 120081 ) who showed that the FR4 filament is a lower-amplitude oscillator than 



the other filament in the December 2006 flare event. 
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The spatial features of the velocity power in the different frequency bands show inter- 
esting dynamics, particularly for the filaments. We observe different sections of the filaments 
fluctuating over various timescales, indicating that either the damping of the oscillations is 
spatial dependent or the excitation is not constant over the entire length of the filaments, or 
some combination of both effects. More detailed analysis of this is underway. 

There is both new and complementary information when comparing the original full- 
disk time series and the FFAMs (the movies associated with Fig. [TJ and the FFAMs of 
Fig. IHJ). For example, for the 2005 event, in the raw intensity time series one clearly sees 
the reaction of some of the far-away filaments after flare eruption, but the FFAM indicates 
that some filaments are actually not affected by the (unknown) trigger mechanism, and also 
directly provides the oscillation frequencies of those which are. In the 2006 flare the Moreton 
wave is easily visible in both intensity and Doppler data, yet the response of the filaments 
is rather muted when viewed in the FFAM, suggesting that the oscillation activation was 
not as effective. Certainly more events need to be treated in such a systematic fashion to 
understand if these are common occurences or not. 

We finally point out that at the photosphere, active regions typicall y show an enhance - 



ment of acoustic power at higher frequencies, known as acoustic halos (IBraun et al.lll992l ). 
This is quite different than the case here where we see higher amplitudes in the lower- 
frequency bandpasses. How the frequency content of waves changes through the different 
layers of the atmosphere could help in understanding the local environment parameters that 
are quite difficult to directly observe, such as magnetic field and density. EUV waves can 
also presumably interact with filaments, although for the two events studied here there were 
no EUV data that could be used to identify any such contribution. In the future this method 
will be used to study a range of various chromospheric data. 



We acknowledge fruitful conversations with Stuart Jefferies, Michael Kirk, Frank Hill, 
and R.T. James McAteer. This work was supported in part by an NSF PAARE award 
AST-0849986. 
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Fig. 1. — Snapshots (pre-flare) in Ha intensity from the ISOON telescope for the (a) 2005 
M-class flare and the (b) 2006 X-class flare. Three regions of interest are marked for each 
event, where "FR" denotes a flaring region and "FL" denotes a filament. See online version 
for movies of the intensity and Doppler time series. You can view all movies online at 



http: / / astronomy.nmsu.edu/jasonj/Halpha_movies/ 
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Fig. 2. — The two flaring regions. Each panel is a snapshot of an FFAM, showing the 
logarithm of the amplitude computed from 60 min of Ha velocity data for the 1, 3, 5, and 
7 niHz frequency bands, from left to right. The top (bottom) row shows the amplitude in 
the flare region FRi (FR2) from segments centered at about 9 and 23 minutes after flare 
initiation, respectively. The color scale is clipped at 0.5 of the maximum velocity amplitude 
to emphasize smaller scale features. See online version for full FFAMs of both intensity and 
velocity for each frequency and flare event. 
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Fig. 3. — Filaments during the May 13, 2005 flare. Each panel is a snapshot of an FFAM, 
showing the logarithm of the amplitude computed from 60 min of Ha velocity data for the 
1, 3, 5, and 7 mHz frequency bands, from left to right. The top (bottom) row shows the 
amplitude in two quiescent filament regions FLi (FL 2 ) for each frequency passband, about 
55 and 62 minutes after flare eruption, respectively. The color scale is clipped at 0.5 of the 
maximum velocity amplitude. See online version for full FFAMs. 
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Fig. 4. — Filaments during the December 6, 2006 flare. Each panel is a snapshot of an 
FFAM, showing the logarithm of the amplitude computed from 60 min of Ha velocity data 
for the 1, 3, 5, and 7 mHz frequency bands, from left to right. The top (bottom) row shows 
the amplitude in two quiescent filament regions FL 3 (FL 4 ) for each frequency passband, 
about 35 minutes after flare eruption. The color scale is clipped at 0.5 of the maximum 
velocity amplitude. See online version for full FFAMs. 
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Fig. 5. — Analysis of flare FR^ using standard methods, (a) A snapshot of intensity around 
the flare. The red circle indicates the pixels (radius of 10 pixels) that were averaged to 
compute the time series in panel (b). The horizontal and vertical lines denote slices used for 
the plots in Fig. [6j (b) Time series of normalized intensity and Doppler velocity for the set 
of pixels in the flare region as a function of time, (c) Wavelet power spectra of the velocity 
and intensity light curves from (b). The (linear) gray scale in these plots is such that black 
denotes higher signal, and has been clipped to half of the maximum power to emphasize 
details. 
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Fig. 6. — Stack plots of flare FRi along the slices indicated in Fig. [5k- The top two panels 
show the intensity and velocity time dependence along the A-B vertical line through the 
flare. The bottom two panels show the intensity and velocity time dependence along the 
C-D horizontal line through the flare. 
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Fig. 7. — Analysis of filament FLi using standard methods, (a) A snapshot of intensity 
around the filament. The red circle indicates the pixels (radius of 10 pixels) that were 
averaged to compute the time series in panel (b). The horizontal and vertical lines denote 
slices used for the plots in Fig. [HJ (b) Time series of normalized intensity and Doppler 
velocity for the set of pixels in the filament region as a function of time, (c) Wavelet power 
spectra of the velocity and intensity light curves from (b). The (linear) gray scale in these 
plots is such that black denotes higher signal, and has been clipped to half of the maximum 
power to emphasize details. 
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Fig. 8. — Stack plots of filament FLi along the slices indicated in Fig. [7k. The top two 
panels show the intensity and velocity time dependence along the A-B vertical line through 
the filament. The bottom two panels show the intensity and velocity time dependence along 
the C-D horizontal line through the filament. 
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Fig. 9. — Snapshots of full-disk velocity FFAMs around peak flare intensity. Panels (a) and 
(b) are the 1 mHz and 3 mHz bands for the 2005 flare, respectively. Panels (c) and (d) are 
the 1 mHz and 3 mHz bands for the 2006 flare, respectively. All panels are plotted on the 
same gray scale. See the full-disk FFAMs online. 



